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C<1 : ABSTRACT 

Q> . In this lecture I will review some recent progress in improving the accuracy of 

I the calculation of density perturbations resulting from inflation. 

1. Introduction 

The early universe was very nearly uniform. However, the important caveat in that 
statement is the word "nearly." Our current understanding of the origin of structure 
^ I in the universe is that it originated from small "seed" perturbations, which over time 
■ grew to become all of the structure we observe. The best guess for the origin of these 
perturbations is quantum fluctuations during an inflationary era in the early universe. 

The basic idea of inflation is that there was an epoch early in the history of the 
universe when potential, or vacuum, energy dominated other forms of energy density 
such as matter or radiation. During the vacuum-dominated era the scale factor grew 
exponentially (or nearly exponentially) in time. In this phase (known as the de Sitter 
phase), a small, smooth spatial region of size less than the Hubble radius at that time 
can grow so large as to easily encompass the comoving volume of the entire presently 
observable universe. 

If the early universe underwent this period of rapid expansion, then one can un- 
derstand why the universe is approximately smooth on the largest scales, but has 
structure (people, planets, stars, galaxies, clusters of galaxies, superclusters, etc.). 
Inflation also predicts that the cosmic background radiation should be very nearly 
isotropic, with small variations in the temperature. Perhaps all of the structure we see 
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in the universe is a result of quantum-mechanical fluctuations during the inflationary 
epoch. In this lecture I will explore this possibility. 

Because nearly all of the students are familiar with the basics of cosmology, I will 
not bother to define familiar terms and notation. In general, the notation follows that 
in The Early Universe (Kolb and Turner, 1990), except that here the scale factor is 
denoted by a{t). 

Since this is a school, I will not provide an exhaustive list of references to original 
material, but refer to several basic papers (including several review papers) where 
students can find the references to the original material. The list of references include 
Bardeen (1980); Stewart (1990); Mukhanov, Feldman, and Brandenberger (1992); 
Liddle and Lyth (1993); and Lidsey, Liddle, Kolb, Copeland, Barreiro, and Abney 
(L^KCBA) (1997). 

2. Evolution of Perturbations 
2.1. Life Beyond the Hubble Radius 

An important part of this lecture will be the interplay of physical length scales with 
the Hubble radius. The time-dependent Hubble radius is defined as the inverse of 
the expansion rate: Ruit) = H~^(t) = [87rGp(t)/3]~^/^ (the last part of the equation 
comes from the Friedmann equation for a spatially fiat universe). In a radiation- 
dominated (RD) universe p oc and in a matter-dominated (MD) universe p oc a~^, 
so Rh oc in an RD universe and Rh oc a^^^ in a MD universe. 

First, let us review what is meant by "crossing" the Hubble radius. For the sake 
of illustration, let's take a length scale A to be at present Aq = 300/i~^Mpc. Today 
the Hubble radius is Rnito) = Hq^ ~ 3000/i-iMpc, so Xo/Rnito) = 10"^ and Aq 
is said to be "within" the Hubble radius today. Any physical length scale increases 
in proportion to the scale factor in an expanding universe. The scale which is Aq 
today, was smaller in the early universe by a factor of a{t)/ao = 1/(1 + z), where Qq 
is the present scale factor. The Hubble radius also depends upon a{t), e.g., Rh = 
[a(t)/ao]^/2 = l/{l + zf/'^ in the MD era. So during the MD era the ratio X{t)/RH{t) 
depends upon redshift as X{t)/RHit) = [Ao/i?//(to)](l + z)^^^ = I0'\l + zf''^- So 
for z > 100, the length scale A was outside the Hubble radius, for z < 100, the 
length scale A was inside the Hubble radius. At 2; = 100 we say that a length scale 
of 300/i~^Mpc crossed the Hubble radius. 

Since X{t) / Ruit) decreases in time in a radiation-dominated or matter-dominated 
universe, any physical length scale A starts larger than Rh-, then crosses the Hubble 
radius (A = H^^) only once. This behavior is illustrated by the left side of Fig. 0. 

In order for structure formation to occur via gravitational instability, there must 
have been small preexisting fluctuations on physical length scales when they crossed 
the Hubble radius in the RD an MD eras. In the standard big-bang model these 
small perturbations have to be put in by hand, because it is impossible to produce 
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Figure 1: The behavior of the Hubble radius, Rh, and the physical size, A, with the 
scale factor a during normal expansion (left) and inflationary expansion (right). 



fluctuations on any length scale while it is larger than Rh- Since the goal of cosmology 
is to understand the universe on the basis of physical laws, this appeal to initial 
conditions is unsatisfactory. 

That any length scale crosses Rh only once is not a fundamental result of anything 
sacred like Einstein's equations, the cosmological principle, or special relativity, but 
it depends upon the assumption of the equation of state. To see how changing the 
equation of state changes the ratio X(t) / Rnit), let's define L(t) to be the dimensionless 
ratio \{t)/RH{t). Obviously, if L{t) is smaller than unity, the scale is within the 
Hubble radius and it is possible to imagine some microphysical process establishing 
perturbations on that scale, while if L{t) is larger than unity, it is beyond the Hubble 
radius and no microphysical process can account for perturbations on that scale. Now 
Rnit) = H-^it) = a{t)/a{t) and X{t) oc a{t), so L{t) is proportional to a(t), and L{t) 
scales as a(t), which from the Einstein equation is proportional to — (p + 3p). There 
are two possible scenarios for L(t) depending upon the sign of p + 3p: 
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Rnit) grows faster than X(t), 



occurs for p + 3p > 



Rnit) grows more slowly than \(t), occurs for p + 3p < 0. 



If during some epoch the equation of state was such that p + 3p < 0, then scales 
larger than Rh remained larger than Rh, while scales smaller than the Hubble ra- 
dius were destined eventually to grow larger than the Hubble radius. The opposite 
behavior obtains during the standard RD and MD epochs when p + 3p > 0. During 
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Figure 2: In the reference unperturbed universe, constant-time surfaces have constant 
spatial curvature (zero for a flat FRW model). In the actual perturbed universe, 
constant-time surfaces have spatially varying spatial curvature. 

these epochs scales smaller than Rh remain smaller than Rh and scales larger than 
Rh eventually become smaller than Rh- 

Now if p + 3p < in the early universe and p + 3p > in the later universe, then 
it is possible to have a "double-cross" situation illustrated on the right side of Fig. 
|I|. In the double-cross scenario, length scales start smaller than the Hubble radius 
during the phase when p + 3p < (the inflationary phase), cross the Hubble radius, 
then remain larger than the Hubble radius. During the standard phase, scales of 
astrophysical interest start larger than the Hubble radius, cross the Hubble radius, 
then remain smaller than the Hubble radius. 

Unlike the standard model, the double-cross model has the feature that it is 
possible to imprint perturbations on a scale as it crosses the Hubble radius during 
the inflationary phase, so one can imagine a reason to have preexisting perturbations 
on scales recrossing the Hubble radius during the RD-MD epochs. 



4 



2.2. Metric Perturbations on Scales Larger than Rh 



What we are interested in following the evolution of a spacetime which is neither ho- 
mogeneous nor isotropic. We will do this by following the evolution of the differences 
between the actual spacetime and a well understood reference spacetime. So we will 
consider small perturbations away from the homogeneous, isotropic spacetime (see 
Fig. 0). 

When one studies "perturbations" it is necessary to specify a reference background 
system. The reference system in our case is the spatially flat Friedmann-Robertson- 
Walker (FRW) spacetime, with line element ds"^ = a^(r) {dr"^ — Sijdx'^dx^}, where r is 
conformal time, related to "normal" time by a'^dr'^ = dt^. Sometimes equations will 
be written in terms of conformal time r, and sometimes in terms of coordinate time 
t. Derivatives with respect to conformal time will be denoted by a prime, while usual 
time derivatives are denoted by a dot, e.g., the Hubble parameter can be defined as 
H{t) = d/a, or ?i(T) = a' /a = Ha. 

The most general form of a metric describing small perturbations away from the 
fiat FRW metric contains scalar, vector, and tensor perturbations [the covariant de- 
composition of Sgfj_iy is given in Stewart (1990)]. For the moment, we will only be 
interested in the scalar perturbations. The perturbed line element including the 
scalar perturbations can be written in terms of four scalar functions {A, B, ip, E}: 

ds^ = a2(r) {(1 + 2A)rfr2 - 2diB dx'dr - [(1 - 2^lj)6,j + 2didjE] dx'dx^} . (2) 

Now because of the residual gauge freedom, not all of the four scalar perturbation 
functions {A, B, ip, £"} are independent. For instance if one works in the synchronous 
gauge, all hypersurfaces have the same time. In this gauge A = B = 0, and the line el- 
ement is c/s^ = ^^(t) {dr^ — [(1 — 2ip)5ij + 2didjE] dx^dx^}. In the longitudinal gauge 
B = E = 0, and the line element is ds^ = a'^^r) {(1 + 2A)dT'^ - [(1 - 2i!)Sij] dx'dx^}. 

It is sometimes bewildering to read the literature because everyone seem to have 
his/her favorite gauge. But really smart people support freedom of choice, and work 
with combinations of the gauge-invariant scalar functions and $ first found by 
Bardeen (1980): 

^ = tjj-n{B-E') 

$ = A + a-^[{B - E')a]'. (3) 
Note that in the longitudinal gauge ^ = A and = 

2.3. Perturbations in the Stress-Energy Tensor 

Inflation assumes that the universe was dominated by something with an equation of 
state that satisfies the inequality p + 3p < 0. Such a component of the energy density 
is usually called "vacuum energy." Since today we know that the vacuum energy 
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must be very small (compared to what is required for inflation) ,[] any inflationary 
model has to have some dynamics for changing the vacuum energy. It is convenient 
to imagine that the dynamics of the change in the equation of state during inflation 
is described by the usual dynamics of a minimally coupled scalar fleld evolving under 
the influence of a scalar potential. This mysterious scalar fleld, denoted by (j), is 
known as the inflaton, and its potential, V{(f)), is known as the inflaton potential. 

One assumes that the inflaton fleld is homogeneous in the reference spacetime, 
0(x, r) = 0o('7")5 and satisfles the equation of motion (pQ + ?>H(j)Q + V,,^ = 0. (Since 
"prime" was used to denote d/dr, I will use X,^ to denote dX/dcj).) This fleld equa- 
tion, together with the Friedmann equation, can be solved to flnd the evolution of 
the background spacetime and scalar fleld. Alternatively, one can view the scalar 
fleld itself as the dynamical variable of the system. This allows the Einstein-scalar- 
fleld equations to be written as a set of flrst-order, non-linear differential equations 
(Grishchuk & Sidorav 1988; Muslimov 1990; Salopek & Bond 1990, 1991; Lidsey 
1991a) 



m 



in 



H,^. (4) 



In the actual perturbed spacetime there are small perturbations about the back- 
ground value: 0(x, r) = 0o(t) + 50(x, r). Of course we will be interested in the 
evolution of 6(f). 

Now just as the metric perturbations are gauge dependent, the scalar fleld fluc- 
tuations are also. But one can circumvent the usual problems associated with gauge 
freedom by constructing a suitable gauge-invariant scalar fleld fluctuation, 50 = 
50 + 0'o(i?-E'). 

2.4. Perturbation Spectra 

Of great convenience is the particular gauge-invariant quantity 

7^ = ^ - ^50 = V - 5^50. (5) 

Clearly 71 deflned by the flrst equality of Eq. (^) is gauge invariant because it is 
constructed explicitly from gauge-invariant terms. However even the second form in 
Eq. (H) is gauge invariant, as when one performs a gauge transformation the non- 
gauge-invariant terms in i/j cancels the non-gauge-invariant terms in the 50 term. 

TZ has a simple physical interpretation in the synchronous gauge, where V^TZ = 
^R/A with the three-dimensional Ricci curvature on the spatial hypersurface. Now 

^And tasteful people assume that today the vacuum energy vanishes, i.e., A = 0. 
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Figure 3: Quantum fluctuations in the infiaton field during inflation lead to fluctua- 
tions in the gauge-invariant quantity TZ. As a particular length scale passes outside 
of the Hubble radius, the fluctuations are "frozen in" because 71 is constant outside 
the Hubble radius. When the length scale reenters the Hubble radius in the RD or 
MD era, fluctuations in TZ appear as energy density fluctuations. 

the usefulness of TZ follows from the fact that as shown by Bardeen (1980), 71 is 
constant on scales much larger than Rh- 

The picture of the generation of quantum fluctuations during inflation can be 
appreciated by studying Fig. ^. 

Now 7Z is related to the observationally determined power spectrum. The first 
step in developing the relation is to expand 7Z in terms of Fourier modes 7lk 



the power in 71 per decade of k. If the curvature perturbation is independent of k, 
then the "power-per-decade" is constant, and \7Zk\ oc k'^/"^. Putting in the factors of 
27r, we define the scalar spectrum As{k) by0 



where As{k) is the primordial scalar density perturbation power spectrum. If 7Z is 
independent of scale outside of the Hubble radius, then As{k) will be independent 

^The exact constant of proportionality is a matter of convention, see L^KCBA. 




(6) 





(7) 
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of k. The primordial power spectrum, As{k), is related to Ps{k), the power spec- 
trum observed in large-scale structure (LSS) surveys and cosmic background radiation 
(CBR) experiments. 

To find the relation between As{k) and Ps{k), it is important to appreciate that 
As{k) is the amplitude when a scale k crosses the Hubble radius, i.e., when k = aH. 
Now if we specify the perturbation spectrum on a particular space-like hypersurface, 
rather than as each scale crosses the Hubble radius, we have to realize that we are 
specifying a gauge-dependent quantity beyond the Hubble radius. We will denote the 
perturbation defined this way as {6p/p)k- In the synchronous gauge and in the co- 
moving gauge, the density perturbation of wavenumber k grows as {Sp/p)k oc {aH)~^ 
for k < aH in both the MD and RD eras. So for scales well outside the Hubble 
radius, {6p/p)k oc [k'^ / {aHy]As{k) , so that when aH = k, {6p/p)k = As{k). 

For scales inside the Hubble radius the synchronous gauge and the comoving 
gauge coincide, and {Sp/p)k is approximately constant in the RD era and grows as 
{6p/p)k oc {aH)~^ in the MD era. So just around the time of matter domination, 
on scales smaller than Rh (i.e., k ^ {aH)EQ), {6p/p)k has the approximate value it 
had when it crossed the Hubble radius, so {6p/p)k ~ As{k) for k ^ {aH)EQ- After 
matter domination {5p/p)k grows as {aH)~'^ ~ a oc t"^^'^ on all scales, so {5p/p)k will 
continue to have the shape it did just after matter domination (at least in the regime 
of linear evolution) . 

The transition between scales larger than Rh at tsQ and scales smaller than 
Rh at tEQ can be encoded in a "transfer function" T{k), by writing {6p/p)k = 
{k/ aHYT{k)As{k) (see e.g., Liddle and Lyth, 1993). In order to reproduce the 
behavior discussed above, the function T{k) must have the limiting forms T{k) 1 
for k <t:aH and T{k) fc"^ for k > aH. 

Now the power spectrum Ps{k) is defined by {5p/ p)\ oc k^Ps{k), so in terms of 
the primordial spectrum As{k) and the transfer function T{k), Ps{k) is given by 
Ps{k) oc kT'^{k)A'g{k). Note that if the primordial spectrum is independent of scale, 
i.e., if A'g{k) is independent of k — the Harrison-Zel'dovich spectrum, then Ps{k) oc k 
for k <^ {aH)EQ and Ps{k) oc k~^ for k ^ {aH)EQ- If we write As{k) as a power 
law, Ag{k) = 74|(A;o)(fc/fco)'^~^) then Ps will be a power law also: Ps{k) oc k"', with 
n = 1 corresponding to the value for constant amplitude perturbations at Hubble 
radius crossing. 

Finally, we must understand the relation between wavenumber k and field value 
(J). During the evolution of the scalar field the background value of changes in 
time. Now associated with a particular value of is a length scale with comoving 
wavenumber k crossing the Hubble radius at the time the scalar field value is (p. The 
easiest relation to find is the differential form found from the expression k = aH: 



dink _ H,tf) a,tf) _ if,^ 



An H 



8 



d(j) H a H 




where the last equality follows from Eq. (^. 
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3. Scalar Perturbations from Inflation 



3.1. Textbook Treatment 

In the standard textbook treatment one expands fluctuations of the inflaton field in 
a Fourier expansion 

{6<pY' = /f ^'I^^I'- (9) 

Then one identifies {5(f))k = k^^'^\5(f)k\ as the fiuctuation in the infiaton on the length 
scale 27r/k. One knows that a scalar field has quantum fiuctuations in deSitter space, 
or more precisely, the quantum fluctuations of a scalar field in deSitter space differ 
from the quantum fiuctuations of a scalar field in fiat space. The quantum fiuctuations 
of a field in deSitter space at Hubble radius crossing (i.e., on scales k = aH, where k 
is the comoving wavenumber) result in {6(f)) k=aH = H/27i. 

Now in the ip = ^ gauge, using the above result for the scalar field fiuctuation 
gives 71 = {H/(j)o){S(j))k=aH ~ H^/<Po- Now we can express and 0o in terms of 
the infiaton potential and its derivative. The background equation of motion from 
Section 2.3 in the slow-roll limit (ignoring the term) gives (p ~ V,^/3H, and the 
Friedmann equation relates H and V{(f)): = Svrp/Smp; ~ \^(0)/mp;. Substituting 
V and V,(^ gives the familiar result for the perturbation spectrum first found in this 
manner by Bardeen, Steinhart, and Turner (1983), 

e/"JZ, ~ As{k) ~ (10) 

Since the scale factor increases so rapidly during infiation, all astrophysical scales 
of interest correspond to a rather narrow range of infiaton field values. For fiat 
potentials, V{(f)) and V,(j){(f)) does not change much during infiation, so one expects 
As{k) to be roughly independent of k. This is the reason for the often repeated 
"result" that infiation leads to an approximate Harrison-Zel'dovich spectrum of scalar 
density perturbations. 

3.2. The Three Step Program for Better Predictions 

The calculation of As{k) in Eq. (p!OD was sufficiently accurate for a decade. But even 
with present-day data, and especially looking forward to the wealth of information 
expected in the near future, such as the angular power spectrum of CMB fiuctuations 
up to multipole number of more than 10'^, more accurate predictions are required. 

As the result of considerable effort, in the past few years some progress has been 
made in improving the accuracy of the calculation of the density perturbation spec- 
trum. Let me describe three basic steps in the road toward better accuracy: 

1. a better treatment of the background classical dynamics by use of Hamilton- 
Jacobi formalism. 
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2. a better formalism of quantum corrections by use of the variational approach, 
and 



3. the calculation of the spectra in terms of slow- roll parameters. 



I have already discussed the advantage of treating H as fundamental, and pa- 
rameterizing its evolution by (f) rather than time. In principle, the Hamilton- Jacobi 
formalism enables one to treat the dynamical evolution of the scalar field exactly, 
at least at the classical level. In practice, however, the separated Hamilton- Jacobi 
equation, the first line of Eq. (|^), is rather difficult to solve. On the other hand, 
the analysis can proceed straightforwardly once the functional form of the expansion 
parameter has been determined. This suggests that one should view H^cf)) as 

the fundamental quantity in the analysis (Lidsey 1991b, 1993). This is in contrast 
to the more traditional approaches to inflationary cosmology, whereby the particle 
physics sector of the model — as defined by the specific form of the infiaton potential 
V{(t)) — is regarded as the input parameter. 

It proves convenient to express the scalar and tensor perturbation spectra in terms 
of H{(f)) and its derivatives. The slow- roll approximation is an expansion in terms of 
quantities derived from appropriate derivatives of the Hubble expansion parameter. 
Since at a given point each derivative is independent, there are in general an infinite 
number of these terms, but only the first few enter into any expressions of interest. 
The first three are defined as 



m = — 



v{4>) 



m 



PI 



H(f) 



_0_ ^ "^pi H,^^{(f)) 
4n H{(t)) 

1/2 



mpie. 



Vie 



vre 



(0)' 



47r 



m 



1/2 



erj 



pr 



47r 



1/2 



(12) 



(13) 



One need not be concerned as to the sign of the square root in the definition of ^; it 
turns out that only and not ^ itself, will appear in our formulae. We emphasize 
> implies that = 



that the choice > implies that y/e = — ^ m^^/ An H H . 

One can show that infiation ends when e = 1. The slow- roll approximation, as 
I use it here, involves assuming {e, 77, ^} are all less than unity. This is somewhat 
more restrictive than just saying that H changes slowly enough for infiation to occur: 
that only requires e < 1. 

Probably the most important advance is the development of the Mukhanov for- 
malism for the perturbation calculation. Recall that the action for the Einstein-scalar 
field system is 

1 



S 



d Xy 



-9 
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(14) 
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with Qf^y = g^^^ + 5gf,u{A, B, E) and = 0o(^) + '50(x, r). 

Before quantizing the system, one must express the theory in terms of the ADM 
variables, expand to second order in the perturbations, apply the background field 
equations, and integrate by parts when judicious. Now the procedure is quite long 
and tedious. Details can be found in the review article of Mukhanov, Feldman, 
and Brandenberger (1992). When the dust settles, the variation of the action can 
be expressed in terms of the dynamical variable u = a{6(j) + (pQip/Ti.) = zTZ, where 
z = a(l)'^/T-i = acpo/H: 

S2S = ^J d^x (^'^ - 5'^u,i . (15) 

Now this is really remarkable, because the complicated dynamics of scalar field per- 
turbations coupled to metric perturbations can be cast into the dynamics of a system 
we know well: a scalar field u in fiat spacetime with (time-dependent and negative) 
mass m? = —z"/z. 

Now scalar field theory in flat space is well understood. So we can use the tool of 
scalar quantum field theory as a sort of hammer to pound out the answer. Of course, 
we have to make sure we have the right tool, because as the saying goes, "when you 
have a hammer in your hand, everything you see looks like a nail." 



3.3. Quantization 

The quantization of the action in Eq. ([T5|) is really rather straightforward: From the 
scalar field m(x, r), form the conjugate momentum 7r(x, r), and form the Hamilto- 
nian from m(x, r) and 7r(x, r). Then promote the classical field and its conjugate 
momentum to operators, m(x, r) — > {i(x, r), and 7r(x, r) 7f(x, r) and impose the 
canonical equal-time commutation relations [m(x, r), M(y, r)] = [vf (x, r), vf (y, r)] = 0, 
and [S(x, r), 7f(y, r)] = i5^{x - y). 

If one expands the field operator in Fourier modes associated with creation and 
annihilation operators al and dk, then the field becomes 

u(x, r) = (27r)-3/2 J d\ [ufc(r)afce*'^-^ + <(r)4e"*'' ''] . (16) 

The field equation for is the familiar Klein-Gordan equation 

< +{e--\uu = Q. (17) 



Of course we must specify the boundary condition. In our case, we want Uk^aH ^ 
^-ikT I ^ i.e., plane- wave solutions. 

It is pleasing to note that any solution to Eq. ( [T7| ) will have the feature that well 
beyond the Hubble radius, IZk will be constant. Note that in the limit — > the 
field equation becomes u'l, — [z" /z)uk = 0, which obviously has solution Uk oc z. Now 
since TZk = Uk/z, on scales much larger than the Hubble radius TZk constant. 
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4. Tensor Perturbations 



In addition to the scalar perturbations in Eq. (|]), the most general metric contains 
perturbations that transform like a tensor on the spatial hypersurfaces. These tensor 
perturbations enter the metric as 

ds^ = a?{T) dr^ - {5ij + hij)dx^dx^ . (18) 

As can be seen by explicit calculation from the Einstein equations, the metric pertur- 
bation hij does not couple to the stress-energy tensor, but describes the propagation 
of gravitational waves. The gravitational waves are not important for large-scale 
structure, but they do have an effect of the CMB, at least for small multipole num- 
ber. 

Since by construction hij is a transverse, traceless tensor, it has two degrees of 
freedom, usually denoted as and h+. (From the quantum view, gravitational waves 
are the propagating part of the gravitational degrees of freedom, corresponding to a 
massless spin- two particle, which of course has two degrees of freedom.) 

Now just as was done for the scalar degrees of freedom, one substitutes the metric 
Eq. (p!8| ) into the Einstein-Hilbert action, and expands to quadratic order in h^j, with 
the result 

2 

3^3="^ I drdha'iT) d.h'.df'hi'. (19) 
647r J ■' 

Since our goal is quantization, and we know how to quantize scalar field theory, we 

want to make 828 look as much as possible like the action for a scalar field. To this end, 

it is very convenient to define the rescaled variable Pj(x) = (mpj/327r)"^/^a(r)/i*^(x). 

In terms of P^j{x), 62S becomes 



52S=- 



^-Jdr d\ (drP, {drP'^) - 5"^" {dmP, ') {dnP)) + ^P, ' P) 



(20) 



This may be interpreted as the action for two scalar fields in Minkowski spacetime, 
each with effective mass squared —arr/a- We can now proceed with quantization 
exactly as in the scalar case. 

Again, just as in the scalar case for TZ, we perform a Fourier decomposition 
of But since there are two degrees of freedom, we must include a polariza- 

tion tensor e*j(k; A), which satisfies the conditions e^j = eji, en = 0, fc*ejj = and 
e*j(k, A)ej-'*(k, A') = 6xx'- The analysis is further simplified if we choose ejj(— k. A) = 
e*j(k, A). The Fourier decomposition can be written as 

P) = j:Jj^2 v,,xe^^{k-X)e^^-. (21) 

In terms of v^, the spectrum of gravitational waves, Axik), is defined as 



k 



.3 



, 2 M^x) = ^-^-^ ^rik) 5^'\k - 1) . (22) 
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Now returning to the quantization of the perturbations, in momentum space the 
tensor perturbation action is 



60S 



- E 

A=+,x 



drdk 



{dr WAY 



k'^ \vk,x\ 



(23) 



We can now quantize Vk^x in the usual way, promoting the field to an operator with 
canonical quantization conditions. 
The mode equation for Vk becomes 



dr^ 



+ e 



a 
a 



Vk = 0^ 



(24) 



This equation can be compared to the mode equation for scalar perturbations, Eq. 
(p!7|). The mode equation is somewhat simpler than the scalar case because a" /a = 
2a^if^(l — e/2) is generally a simpler function than z"/z. 



5. A Variety of Models (Some Realistic, Others 
Illustrative) 

5.1. Solution Procedure 

The procedure is simple (in principle): solve Eq. (|T7|) for Uk, then find TZk = Uk/z to 
give As{k), which together with a transfer function, yields the power spectrum which 
can be compared to observations. Then solve Eq. (|2^) for Vk, to give AT{k). 

The trouble is that exact solutions to the wave equations are hard to find, partly 
because the mass terms are so complicated: 



- = 2a''H' 

z 

a" , , / 1 

— = 2a'H' 1 - -e 
a V 2 



l + e-^r] + e'~2er] + ^r]' + 



(25) 



where the r dependence of H, e, r], and ^ are found from their dependence upon 
(f). In fact, only two exact solutions of Eqs. (|1^ and (p^) are known. The first is a 
power-law solution found by Stewart and Lyth (1993), and the second, yet unnamed, 
has been found by Easther (1996). 

The first step is to express the conformal time, dr = dt/a{t) in terms of aH and 
the slow-roll parameters. In general the result is 



T 



da 



aH 



+ 



da. 



(26) 
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If e is constant, then r^^ = —aH{l — e) (r is negative during inflation, with r = 
corresponding to the infinite future). If e is not constant, then integrating by parts 
an infinite number of times, one can obtain 

r = + expansion in slow-roll parameters ( etc. , (27) 

aH 1 — e aH 

where ( = e — rj, and e can now have arbitrary time dependence. 

In the next section I will review the exact power-law solution, and and the section 
after that I will discuss how to use that exact solution to construct perturbative 
solutions for other models. 



5.2. Power-Law Inflation 



In the power-law model the Hubble parameter is expressed in terms of the Planck mass 



and a parameter p: H{(j)) oc exp AiKp"^ / pm'pi, which results from a scalar potential of 
the form V{(f)) oc e'^ (Abbott and Wise 1984, Lucchin and Matarrese 1985). Obviously 
this type of potential is not a fundamental, renormalizable scalar potential, but it is 
the type of effective low-energy potential for dilaton-like degrees of freedom in string 
theories and Kaluza-Klein theories. 

For H{(j)) oc e'^, e, i], and ^ will be equal and constant: e = rj = ^ = p~^. 

Now one can proceed to find z"/z and a" /a, with the result z"/z = {v"^ — 1/4) /r^ 
and a" /a = (/i^ - 1/A)/t^, where u = (3/2) + {p - 1)"^ and /i = (3/2) + {p - 1)-^ 
(For power-law infiation u and fi coincide, though in general they do not.) 

For power-law infiation the mode equations are simply a Bessel equation: 



rf2 „ zy2 _ 4-1 
+ ^ 



/i2 - 4-1' 



Mfe = 

Vk = 0, 



(2^ 



which for the boundary conditions we impose are solved by Hl^\-kr) and Hj^^-kr), 
Hankel functions of the first kind of order u and /i. 

We are interested in the asymptotic forms of Uk/z and Vk for k <^ aH, which are 
easily found to be 



Uk 

Vk 



gi(i/-l/2)7r/22!^-3/2 

above with u 



r(3/2) 



i-kr) 



(29) 



which yields As{k) oc H'^/\H'\ and A^ik) oc H, with both expressions evaluated at 
k = aH. Now using the fact that at Hubble radius crossing H{(f)) oc k^/P from Eq. 
(1), we find a power-law spectrum As{k) and A-rik) proportional to k'^/^. 
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The scalar spectral index is defined as ri(A;) — 1 = dlnA^/dlnk. Writing A'g{k) oc 
the above power-law spectrum gives n — I = — 2/p, a departure from the n = 1 
Harrison-Zel'dovich result. Defining the tensor spectral index, nT{k) as nT{k) = 
din A^/ dink, for power- law inflation ut = — 2/p. 



5.3. General Potentials 

After working hard to find an exact solution, we can now make an expansion about 
it for general potentials. The power-law inflation case corresponded to the slow-roll 
parameters being equal, and hence exactly constant. In general they can be different, 
which means they will pick up a time dependence. 

Assuming that e, as well a.s ( = e — r] are small, then Eq. (^) can be approximated 
to give T = -{1 + e)/{aH). 

Having this expression for r, we can now immediately use Eq. (pSf), which must 
also be truncated to first-order. This gives the same Bessel equation Eq. (|^), but now 
with ly given byz/ = 3/2 + 2e — r/ and fi given by/i = 3/2 + e. The assumption that 
treats e as constant also allows r] to be taken as constant, but crucially, e and r] need 
no longer be the same since we are consistent to first-order in their difference. The 
differences between further slow-roll parameters and e lead to higher order effects, 
and so incorporating e and i] in this manner is applicable to an arbitrary infiaton 
potential to next-order. The same solution Eq. (pQ]) can be used with the new form 
of z/, but for consistency it should be expanded to the same order. This gives the 
final answer, which is true for general inflation potentials to this order, of (Stewart 
& Lyth 1993) 

Asik) = <( [1 - (2C + l)e(0) + Cr^(0)] ^ ^^'^^ 



5v^ [ rripi 



2 i^(0)l 



k=aH 



Arik) = <{ [1 - (C + l)e(0)] ^ ^ . (30) 

k=aH 



by/n \ mpi j 

where C = — 2-(-ln2-|-7 ~ —0.73 is a numerical constant, 7 being the Euler constant 
originating in the expansion of the Gamma function. Of particular interest is the 
ratio 

-^ = e(0) [l + 2C{ei<P)-v{m • (31) 

It is useful once again to point out that the < — > k connection is made through 
Eq. (P), which can be written in the form 



dink I An e-1 ^3^) 



m 



PI 



For the spectral indices n{k) — 1 = In v4|(fc) /dink and nT{k) = din A'^{k) /dink, 
it is easy to show that 

n{k) = l-Ae + 27]-[8{C + l)e^-{Q + 10C)e7] + 2C^'^], 
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observable 


lowest-order 


next-order 


AUko) 




H{<i)o), e(0o) 




H{<Po), e(0o) 


^(0o), e(0o), r/((/'o) 


Al{ko)/Al{k,) 


e(0o) 


e(0o), ^(0o) 


nriko) 


e(</»o) 


e(0o), v{(Po) 


1 - n{ko) 


e(0o), r/(0o) 


e(0o), ?7(0o), C(0o) 



Table 1: The observables, A"^, A|, n, and at the point ko may be expressed in 
terms of H and the slow-roll parameters at the point (f)Q. This table lists the inflation 
parameters required to predict the observable to the indicated order. 



ririk) = -2e[l + (3 + 2C)e- 2(1 + C)r;]. (33) 

Obviously, the usual Harrison-Zel'dovich result n = 1 is obtained if the slow-roll 
parameters {e, 1], ^} are all much less than unity. But recall that e = 1 defines the 
end of inflation, so there is no reason to assume that the slow-roll parameters must 
be much less than unity 50 e-folds from the end of inflation. 



5.4. The Consistency Relation 

Before turning to specific models, it is important to recognize a "consistency" relation. 
The overall amplitude is a free parameter determined by the normalization of the 
expansion rate H during inflation (or equivalently the scalar field potential V). On 
the other hand, the relative amplitude of the two spectra is given to lowest order by 

Thus, to lowest order in the slow-roll parameters, there exists a simple relationship 
between the relative amplitude and the tensor spectral index: 

UT = -2^ . (35) 

This is the lowest-order consistency equation and represents an extremely distinctive 
signature of inflationary models. It is difficult to conceive of such a relation occurring 
via any other mechanism for the generation of the spectra. 

Since it is possible for the spectra to have different indices, the assumption that 
their ratio is fixed can be true only for a limited range of scales, but the correction 
enters at a higher order in the slow-roll parameters. 
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5.5. Other Models 



Here I briefly give some results to lowest order in the slow-roll parameters for the 
spectral index in a couple of well-studied inflation models. I will work out polynomial 

chaotic inflation in detail, and only describe the other models and give the results. 

Now in this section we are treating the potential as input, so it is useful to have 
the lowest-order results for the slow-roll parameters in terms of V. These were studied 
by Kolb and Vadas (1994), with the result 

I will use the lowest-order result n = 1 — 4e + 21] and = —nT/2 = e. 

Now of course the slow-roll parameters are a function of 0, which imphes they are 
a function of k. But since we are working to lowest order, we can assume that the 
spectral indices are constant, and the values associated with Hubble radius crossing 
about 50 e-folds from the end of inflation. Generally we will have to flnd the value of 
the fleld 50 e-folds from the end of inflation. We will denote this as 05o. 

The end of inflation is deflned by e(0) = 1, and the deflnition of the number of 
e-folds from the end of inflation is 



N{(t)) = / — = / H dt^ -^d(t)^ / ^= . (37) 

Ja a Jt J<t> (f) \ rripi J4, Je{(f)) 

5.5..1 Polynomial Chaotic Inflation 

Probably because of its simplicity, the most popular inflation model is polynomial 
chaotic inflation, where the potential is assumed to be V{(f)) = (jf. A potential of this 
form has been championed by Linde. 

With the potential in this form, the first two slow-roll parameters are 

and ^ ^ _!!^ ^ + I^i fcH . (38) 

The end of inflation is found by setting e = 1, which gives 0e„d/'^pi = p^/lGir. For 
this model A^(0) = 47r(02 — (f)lnd) /pm-pi, which gives e = pmp;/167r05o = p/{p + 200), 
and 77= (p-2)/(p + 200). 

Finally, the above values of e and 77 give 

n = 1 - 4e + 2r; = 1 ~ 1 - 

' p + 200 100 

A^/Al = e = ^ — (39) 

This gives a flavor of the calculations that can easily be done for the other inflation 
models discussed in the following subsections. The results are given in Table 1. 
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model 


n 


-/ / 


pU W t^l -idW . 






£yvT\ * /l fl'TTz/i^ / 'rim 
Ca.^ a / ±U/I (ji/ / JJH ipi 




1 I'n 


p — 


fi 


9 


p = 10 


0.8 


0.1 


L/*J1 y liVJliildl L/liClVj tlL^ . 








1 _ _i_ 9V1 nn 


p/ i^p -t- ZUUJ 


Tl — 9 


n Qfi 


n ni 


77 — 4 
P ^ 


Q4 


n 02 

yj . 


natural 






1 ± cos(0//) 


1 - mlJST^p 


(m|,;/167r/2) exp(-15m|,,/27r/2) 


/2 = bmlJ^T, 


0.8 


6 X 10"^ 


P = bmlJlQi: 


0.6 


8 X 10-^2 








{l-exp[-(167r/3m|,,)V2^]}2 


0.96 


10-3 


CDM {V = ????) 


1 






Table 2: Well studied inflation models to lowest order. The result A^/Ag = (1 — n)/2, 
true for power-law inflation, is often (incorrectly) used as a general result. The relative 
contribution to tensor modes to the CMB power spectrum for small multipole number 
is approximately d.bA^/Ag. 



5. 5.. 2 Power Law Inflation 

We have already discussed the power- law inflation model. In that model e — rj — 1/p. 
Of course the fact that e is a constant means that some other machinery must be 
introduced for the highly desirable result of an end to inflation. 

5. 5.. 3 Natural Inflation 

Natural inflation is a local Fermilab favorite (Freese, Frieman, and Olinto, 1990). In 
this model the potential takes the form of the potential for a pseudo-Nambu-Goldstone 
boson: 

1/(0) = [1 ± cos{<j>/f)] , (40) 
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where A and / are mass scales. The mass scale / corresponds to the scale of the 
breaking of the original U{1) symmetry, and A is the mass scale associated with an 
explicit breaking term. It is attractive to consider / to be of order m pi and A of order 
the GUT scale. 

Natural inflation is a great example of a model with a non-renormalizable scalar 
potential. Even though the underlying theory may be renormalizable, there is no 
reason to expect that the effective low-energy inflaton potential should be restricted 
to be of a renormalizable form. 



5. 5.. 4 Inflation 



inflation is actually the first model for inflation (Starobinsky 1980). In this model 
the inflaton potential is not a fundamental scalar field, bt has an origin in the gravity 
sector. If one adds a term quadratic in the Ricci scalar to the Einstein-Hilbert action, 



IGTrG 



d X 



R + 



6M2 



(41) 



then by means of a conformal redefinition of the metric the R? term can be elimi- 
nated. Or, more precisely, the extra degree of freedom can be rewritten to look like 
a minimally coupled scalar field with action 



1 



where 



IGTrG 



V{ct>) 



Jd^x^R+Jd^x^ y-r^'d^cpdA - v(0) 



327r 



{l - exp [-(167r/3m|;)^/^0]} 



(42) 



(43) 



Here is an example of an effective inflaton potential where the scalar fleld need 
not be regarded as a fundamental scalar field degree of freedom. This suggests that 
the scalar field analysis described in this paper may be useful for a class of models 
larger than just scalar field models. 



6. So What's Your Point? 

In this lecture I have tried to make several points: 

1. In one-field, slow-roll models of infiation it is possible to make sufficiently accurate 

predictions of the observable parameters such as ^45, At, N, and ut- 

2. The restriction of "one-field, slow-roll" may not be as restrictive as first imagined, 

because many models of inflation can be written in this way even if they do not 
involve a fundamental scalar field to start with. 
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3. Different models make different predictions for the observables. Soon it will be 

possible to sort through the models and start weeding out those not in agreement 
with observation. 

4. There is a consistency relation for these models, although it may be difficult to 

check observationally. 

5. Although not discussed in this lecture, with a little work one can rework the 

expressions for the observables to express the potential in terms of the observ- 
ables. Therefore, one might be able to glean some information about a scalar 
field potential at energy scales of about lO^^GeV from astronomical observations 
(Copeland et al. 1993a, 1993b, 1994; Lidsey et al. 1997; Turner 1993a, 1993b). 
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